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Equation 
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=  c(x,  z),  instantaneous  concentration  of  dif¬ 

fusing  particles  at  (x,  z^' 
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=  deposition.  [gm/(m‘^sec)] 

(7) 

D 
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2  , 

=  maximum  deposition.  [gm/(m  sec)] 

(11), (21), (27) 

D 
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=  deposition  for  zero  atmospheric  diffusion 

(19),  (28) 
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=  H/t)  =  dominant  length  scale 
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f* 

=  H/(ti  + 

(24),  (38) 
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=  instantaneous  flux  of  particles 

(1) 
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=  height  of  source 
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=  von  Karman's  constant  *^0.4 

(6) 

=  proportionality  constants,  function  of  v  only 

(23),  (25).  (28),  (34) 

K 
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(13) 
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which  neglects  atmospheric  diffusion  to 
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approximation. 

(39) 
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r(  ) 
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=  ratio  of  predicted  by  (9),  which  ne¬ 
glects  dispersion  of  terminal  velocities 
to  location  predicted  by  the  analytic 
approximation. 

=  time 

F  mean  wind  speed  at  height  z 

=  friction  velocity 

=  (surface  stress/density) 

=  r  /  u(h)  dh 
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(2) 

(4) 

(6) 

(55) 


=  instantaneous  vector  velocity  of  air  (1) 

=  vertical  component  of  atmospheric  veloc.ity  (4) 

=  distance  downwind  (4) 

=  point  where  maximum  deposit  occurs  (9).  (20),  (26) 


=  vertical  distance  (6) 

=  roughness  parameter  (6) 

s  parameter  determined  by  z  and  depth  (6) 

of  layer  concerned 


=  stability  parameter 


(51) 


=  gamma  function  (7) 

=  dimensionless  deposition  (7) 


=  maximum  deposition,  [dimensionless)  (10),  (21),  (27) 

A^  =  dimensionless  deposition  for  zero  (19) 

atmospheric  diffusion 


A’*'  =  dimensionless  deposition,  incorporating  (24) 

varying  terminal  velocities 

(  =  coefficient  depending  on  atmospheric  (51) 

stability 
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atmosphere 
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Equation 
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(25) 
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(54) 
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=  ratio  of  predicted  by  (21)  or  (30),  which 

neglects  atmospheric  diffusion  to  that  pre¬ 
dicted  by  (27),  the  analytic  approximation. 

(40) 

A*[A*] 

=  ratio  of  predicted  by  (10),  which  ne  ¬ 

glects  dispersion  of  terminal  velocities  to 
that  predicted  by  the  analytic  approximation. 
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(13) 
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(14) 
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(16) 
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Abstract 


\ 


When  a  cluster  composed  of  particles  of  varying  terminal  velo¬ 
city  is  released  in  the  atmosphere,  at  least  two  important  physical 
processes  contribute  to  the  subsequent  separation  of  these  particles. 

The  first  is  separation  of  the  particles  due  to  small  scale  turbulence; 
the  second  is  separation  due  to  the  varying  terminal  velocities  of  the 
particles.  It  is  shown  that  a  term  in  the  equation  of  turbulent  diffusion, 
normally  neglected,  represents  the  contribution  to  the  separation  due 
to  variation  in  terminal  velocity.  This  term  can  be  transformed  by 
conventional  mixing  length  assumptions  and  the  resulting  differential 
equation  may  be  solved  analytically  for  an  adiabatic  atmosphere  pro¬ 
vided  it  is  assumed  that  the  added  term  is  independent  of  the  x 
coordinate.  This  solution  differs  from  that  obtained  by  Rounds  (or 
Bosanquet)  for  the  case  of  monodisperae  particulate  diffusion  in  an 
adiabatic  atmosphere  only  in  that  the  diffusion  coefficient  is  altered 
by  an  amount  which  depends  on  the  statistical  parameters  of  the  dis¬ 
tribution  of  terminal  velocities. 

The  parameters  are  worked  out  in  detail  for  a  lognormal  dis¬ 
tribution  of  terminal  velocities.  By  reference  to  a  symbolic  inte¬ 
gration  of  the  basic  monodisperse  solution  it  is  shown  that,  in  the 
surface  layer  the  pattern  of  non-dimensional  ground  deposition  de¬ 
pends  only  on  the  logarithm  of  the  ratio  of  the  median  terminal  velocity 
to  the  friction  velocity,  and  on  the  coefficient  of  variation  of  the  termi¬ 
nal  velocity  distribution.  Comparison  of  the  analytic  solution  for  the 
polydisperse  cloud  with  exact  numerical  integration  of  the  monodisperse 
solution  indicates  the  error  in  the  analytic  solution  is  tolerable  for  a 
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wide  range  of  the  basic  parameters  despite  the  neglect  of  the  x 
variation  of  the  increment  in  diffusion  coefficient  due  to  the  distri¬ 
bution  of  terminal  velocities. 

The  analytic  polydisperse  solution  is  then  applied  to  a 
diabatic  atmosphere  (after  Godson)  and  criteria  are  developed  to  indi¬ 
cate  when  neglect  of  atmospheric  diffusion  or  particle  terminal 
velocity  distribution  is  permissible.  For  the  range  of  particle  sizes 
of  major  interest  in  air  pollution  studies,  it  is  shown  that  normal  vari 
ations  in  atmospheric  stability  and  wind  speed  are  sufficient  to  make 
neglect  of  either  factor  invalid  for  a  substantial  number  of  occasions. 
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1.  Introduction 

When  a  cluster  composed  of  particles  of  varying  sizes  and 
terminal  velocities  is  released  in  the  atmosphere,  at  least  two  im¬ 
portant  physical  processes  contribute  to  the  subsequent  separation  of 
these  particles.  The  first  is  simply  separation  of  the  particles  due  to 
eddy  diffusion:  the  second  source  of  separation  is  due  to  the  varying 
terminal  velocities  of  the  particles  comprising  the  diffusing  cluster. 

The  relative  importance  of  each  of  these  processes  obviously  depends 
on  the  statistical  distribution,  in  particular  the  mean  and  standard 
deviation,  of  the  terminal  velocities,  and  the  diffusive  power  of  the 
atmosphere  as  expressed  by  some  characteristic  vertical  diffusion 
velocity. 

It  is  the  purpose  of  the  present  paper  to  consider  the  relative 
importance  of  each  of  these  processes  and  to  develop  general  criteria 
for  permissible  neglect  of  either  one  or  the  other.  It  is  convenient  to 
adopt  a  mathematical  model  of  diffusion  from  a  continuous  infinite  ele¬ 
vated  line  source  as  developed  by  Rounds  (1955).  This  treatment  is  mathe¬ 
matically  valid  for  the  diffusion  of  monodispersed  particles  in  the  surface 
layers  of  an  adiabatic  atmosphere  only.  The  line  source  limitation  is  not 
too  serious  since  the  terminal  velocities  of  the  particles  act  in  the  verti¬ 
cal  plane.  The  restriction  to  the  adiabatic  atmosphere  is  more  serious, 
but  Godson  (1958)  has  developed  an  approximate  method  of  applying  the 
original  Rounds  solution  to  the  diabatic  atmosphere. 

The  method  of  procedure  will  be  to  consider  the  mathematical 
solution  for  a  monodispersed  cloud  and  then  to  integrate  this  solution 
over  a  distribution  of  particle  sizes.  The  integration  can  be  performed 


by  numerical  methods  only,  but  the  value  of  the  integral  is  a  function 
only  of  a  virtual  diffusion  velocity  and  the  median  and  standard  devi¬ 
ation  of  the  particle  terminal  velocities.  A  method  of  incorporating 
the  effect  of  varying  terminal  velocity  directly  into  the  diffusion 
equation  is  then  developed,  and  it  is  shown  that  the  resulting  analyti¬ 
cal  solution  is  a  good  approximation  to  the  numerical  integration. 
Criteria,  based  on  the  analytical  solution,  are  then  developed,  in 
order  to  assess  the  relative  importance  of  the  particle  size  distri¬ 
bution  on  the  deposition  patterns  under  specified  terrain  and  wind 
speed. 

2.  Diffusion  of  monodispersed  particles  in  an  adiabatic  atmosphere 
It  is  of  some  importance  in  what  follows  to  consider  the 
general  equation  of  turbulent  diffusion  for  particles  of  uniform 
terminal  velocity.  The  instantaneous  flux  of  particles  will  then  be 

(1)  F  =  (V  +  «)c 

where  V  is  the  instantaneous  vector  velocity  of  the  air,  w  is  the 
instantaneous  velocity  of  the  particle  with  respect  to  the  air,  and  c 
is  the  instantaneous  concentration  of  the  diffusing  particles.  The 
equation  of  continuity  is  now 

(2)  l^s-DivF. 

Upon  averaging  both  sides  of  (2)  and  employing  the  usual  Reynold's 
averaging  postulates,  one  obtains 

(3)  +  Div(  (V  +  c3)c)  ss  -  Div(V'c' 
where  w'  =  w  -  w  ,  c'  *  c  -  c,  and  V  =  V  -  V  . 


+  w'  c ' )  , 


Rewriting  (3)  to  apply  to  the  steady  state  solution  from  an 
infinite  line  source  oriented  normal  to  the  direction  of  the  mean  wind, 
u,  we  obtain 

(4)  u(z)  a  -  Div(w'  +  w')c'  , 

where  w  is  the  vertical  component  of  atmospheric  velocity. 

As  is  customary,  the  term  Div  u'c'  has  been  neglected  in 
comparison  with  the  term  u  SZ/Bx  .  For  a  monodispersed  particle, 
the  term  c  8n/8z  can  be  important  only  near  the  release  line  since 
this  term  is  non-zero  only  for  the  time  or  vertical  distance  required 
for  the  particles  to  achieve  terminal  velocity.  On  the  scale  of  interest 
here,  this  term  is  negligible.  For  convenience  -w  has  been  replaced 
by  fl,  the  terminal  velocity  of  the  diffusing  particles.  Since  the 
particles  are  heavier  than  air  w  will  always  be  negative,  and  f?  will 
therefore  be  a  positive  quantity.  In  the  term  u'c'  on  the  right  hand 
side  of  (4),  u'  represents  fluctuating  departures  from  particle  termi¬ 
nal  velocity.  For  monodispersed  particles,  such  departures  are  due 
to  inertia  and  reflect  the  inability  of  a  particle  to  adjust  instantaneously 
to  the  fluctuating  velocity  field.  Evaluation  of  this  term  obviously  de¬ 
pends  on  the  previous  history  of  the  particle.  It  is  customary  to  neglect 
u'c'  in  comparison  with  the  other  term  involving  w  but  as  we  shall 
see  later,  this  term  becomes  quite  important  when  there  is  a  large 
variation  of  terminal  velocity  amongst  the  diffusing  particles. 

The  right  hand  side  of  (4)  is  transformed  by  conventional 
mixing  length  assumptions,  i.  e.  ,  c'  =  -f  8c’/8z,  u'  =  u*.  K  =  u*f 
where  u*  is  a  characteristic  fluctuating  velocity  and  i  is  a  character¬ 
istic  scale  of  vertical  velocity  fluctuations,  whence 


4 


(5) 

where  for  the  moment  u'i  has  been  set  equal  to  zero.  In  the  surface 

||l 

layers  of  an  adiabatic  atmosphere  i  a  kz  and  u  is  constant  with 
height  whence 


where  u*  is  the  constant  friction  velocity,  equal  to  the  square  root  of 
the  surface  stress  divided  by  the  density;  h  is  a  reference  height, 
identified  as  the  height  of  the  source  in  what  follown;  and  k  is  von 
Karman's  constant  0.4.  The  approximation  to  the  logarithmic  wind 
profile  was  introduced  by  Calder  (1949),  and  employed  by  Deacon  (1949)i 
Rounds  (1955),  F.  B.  Smith  (1957),  and  others.  The  parameter  a  can 
be  determined  as  a  function  of  Zq.  the  roughness  parameter,  and  the 
depth  of  the  layer  concerned  while  u*  can  be  evaluated  from  the  log¬ 
arithmic  law  provided  u(h)  and  z^  are  known. 

Rounds  succeeded  in  solving  (5)  with  the  coefficients  u  and 
K  given  by  (f)  for  the  concentration  of  diffusing  particles  as  a  function 
of  X  and  z.  According  to  this  diffusion  model  the  flux  through  the 
plane  z  a  0  is  simply 

nc(x,  0)  , 

where  c(x,  0)  is  the  concentration  at  z  a  0.  This  quantity  is  interpreted 
as  the  deposition  on  a  flat  surface  at  z  a  0.  The  deposition  at  z  a  0, 
from  a  line  source  at  height  h,  according  to  Rounds  may  then  be 
written 


(7) 


A  fD  _  P  l^\ 

^  ~ -a- ^  <5?) 


e  * 


2 

where  D  is  the  deposit  [gm/(m  sec)],  Q  is  the  source  strength 
[gm/(m  sec]],  T  is  the  gamma  function,  A  is  the  dimensionless  de- 

£2  £2 
P  =  (1  +"a)ku*  =  n 

h  u(h)  ^  H  .  ^^n(h/z  ) 

(1  +  a)\u’*  (1  +  o)^k^ 

The  quantity  r|  has  the  dimensions  of  a  velocity  and  will  be  referred  to 
hereafter  as  the  diffusion  velocity.  The  dimensionless  parameter  p 
therefore  represents  the  ratio  of  particle  terminal  velocity  to  atmo¬ 
spheric  diffusion  velocity.  The  dominant  length  scale  is  f  which  (by 
(6) )  is  the  height  of  the  source  scaled  by  the  ratio  of  the  vertically 
averaged  wind  speed  (determined  from  the  ground  to  h)  to  the  diffusion 
velocity.  Since  the  diffusion  velocity  is  proportional  to  the  wind  speed, 
f  can  be  expressed  as  a  function  of  wind  speed,  wind  shear  and  surface 
roughness  only.  It  will  be  noted  that  f  is  independent  of  the  particle 
terminal  velocity  and  depends  only  on  source  height  and  atmospheric 
wind  and  diffusion  parameters. 

It  may  be  verified  that  the  integral  of  (7)  with  respect  to  x 
from  0  to  90  iu  the  source  strength  Q  and  hence  all  the  material  re¬ 
leased  is  eventually  deposited  on  the  ground.  In  particular  the  maxi¬ 
mum  deposit  occurs  at 


position,  and 

(8) 


(9) 


h  u  (h) 


m  “  1  +  p  "  (1  +  a)I(l  +  a)ku*+  £2] 


and  clearly  demonstrates  that  the  point  of  maximum  concentration  moves 
closer  to  the  source  as  the  terminal  velocity,  (£2),  or  as  the  diffusion 


6. 


velocity  (i\)  increaaee.  It  is  clear  from  (9)  that  insofar  as  the  location 
of  the  maximum  deposit  is  concerned,  the  diffusion  and  terminal  velo¬ 
cities  enter  as  a  linear  combination.  When  il  the  effect  of  the 

terminal  velocities  of  the  particles  is  negligible.  Similarly,  the  effect 
of  diffusion  on  the  location  of  the  point  of  maximum  concentration  be¬ 
comes  negligible  when  >>  il.  Since  (1  +  a)k  is  of  the  order  of  1/2, 
the  two  effects  are  about  equally  important  when  u*  ^  20. 

The  value  of  the  maximum  deposit,  D  ,  fA  ,  in  dimensionless 

m 

units]  is  given  by 


(10) 


A  = 
m 


fD 

m 

"IT 


r(i 


(H-r 


From  (8)  it  is  apparent  that  the  maximum  non-dimensional  deposit 
depends  only  on  the  ratio  of  the  terminal  velocity  to  the  diffusion  velo¬ 
city.  Restoring  dimensions  to  (10)  we  have 


(11) 


-  (1  +  a)0  ^1 

”5“"  hu(h)  r(l  +  p)  V  e  y 


showing  that  the  value  of  the  maximum  deposit  is  proportional  to  the 
terminal  velocity  and  inversely  proportional  to  the  wind  speed  and 
height  of  release. 


3.  Diffusion  of  polydispersed  clouds 

3. 1  Integration  of  monodisperse  solution 

The  above  discussion  indicates  that  there  are  circum¬ 
stances  when  diffusion  is  relatively  unimportant  compared  to  the  terminal 
velocity  of  the  diffusing  particles  and  vice-versa.  We  now  investigate 
the  situation  when  the  diffusing  cloud  is  composed  of  particles  with  a 
distribution  of  terminal  velocities.  Let  the  distribution  of  mass  as  a 


function  of  particle  ternminal  velocity  be  denoted  by  Q(Q). 
Q(f2)  has  the  property 


The 


J  Q(n)  dO  =  Q  . 

0 

If  it  is  assumed  that  particles  with  different  terminal  velocities  diffuse 
independently  of  each  other,  the  deposition  pattern  would  be 


(12) 
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Since  p  is  a  function  of  Q,  (12)  cannot  be  integrated  analytically 
for  any  reasonable  Q(fl). 

We  consider  the  case  now  where  the  distribution  of  mass 
is  lognormal  in  terminal  velocity.  Since  negative  terminal  velocities 
(-fl)  are  excluded,  the  lognormal  distribution  is  a  reasonable  one 
which  may  be  expected  to  arise  in  practice.  For  this  type  of  dis¬ 
tribution 


(13)  Q(n)  =  — exp{-[fnn  -  p]^/2v^}  , 

where  p  is  the  mean  value  of  the  logarithm  of  the  terminal  velocities, 
and  V  is  the  logarithmic  standard  deviation.  The  arithmetic  average, 
fl,  the  median,  fl,  the  arithmetic  standard  deviation,  cr,  and  the  coeffl 
cient  of  variation  (t/O  are  given  by 


(14) 
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Upon  substituting  (13)  into  (12)  and  with  the  help  of  (8)  chang¬ 


ing  the  variable  of  integration  to  p,  (12)  becomes 
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where 


^  a  p  -  inll  a  In(-) 


ri  =  e»"-^ 


(16a) 

and  inversely, 

(16b) 

In  terms  of  the  non-dimensional  variables  employed  above,  the  required 
integral  will  be  a  function  of  x/f,  0,  and  v  only*  The  parameter  0 
represents  the  logarithm  of  the  ratio  of  the  mass -median  terminal 
velocity  to  the  diffusion  velocity.  A  large  negative  value  of  0  implies 
that  the  diffusion  velocity  is  large  compared  to  the  median  terminal 
velocity  of  the  diffusing  particles  while  a  positive  value  of  0  indicates 
that  the  median  terminal  velocity  is  greater  than  the  diffusion  velocity. 

In  effect,  0  estimates  the  relative  importance  of  terminal  velocity 

and  diffusion  velocity  on  the  integrated  pattern  of  deposition.  As 

far  as  the  non-dimensional  pattern  of  deposition  is  concerned,  atmospheric 

diffusion  enters  only  in  the  value  of  0.  Consequently  if  the  logarithm 

of  the  diffusion  velocity  q  is  very  small  compared  to  the  logarithm  of 

the  median  particle  terminal  velocity,  the  pattern  of  non-dimensional 

deposition  is  almost  independent  of  the  atmospheric  diffusion  velocity. 

The  deposition  pattern  for  this  case  will  then  depend  only  on  17  and  v. 

It  should  be  noted  that  since  f  a  H/q,  the  diffusion  velocity  appears  in 
the  denominator  of  the  left  hand  side  of  (15)  so  that  even  when  q  is 
very  small  compared  to  17,  the  non-dimensional  value  of  the  deposition 
will  be  proportional  to  q.  However^he  shape  of  the  curve  of  deposition 
will  not  depend  on  q. 


The  parameter  v  js  the  only  other  independent  parameter  on 
which  (15)  can  depend^  For.  the  case  <<0,  it  is  clear  that  the 
pattern  of  non-dimensional  deposition  will  depend  only  on  and  v. 
The  parameter  v  represents  the  logarithmic  standard  deviation  of 
the  particle  terminal  velocities.  A  more  meaningful  interpretation 
of  V  emerges  from  (14),  namely  the  coefficient  of  variation,  cr/fl  , 
of- the  terminal  velocity  distribution  is  uniquely  fixed  by  the  value 
of  V  .  It  is  this  latter  interpretation  which  will  be  adopted  hereafter. 
For  very  small  values  of  v,  it  seems  clear  that  (15)  should  reduce 
to  (7)  since  the  terminal  velocity  dispersion  around  the  mean  value 
will  then  be  quite  small.  From  these  considerations  it  appears 
that  under  certain  circumstances,  i.  e.  when  T|<<fi,  atmospheric 
diffusion  is  not  an  important  factor  in  describing  the  pattern  of 
ground  deposition  while  for  other  cases,  i.  e.  when  v  or  (r/n 
is  small,  the  spread  of  terminal  velocity  will  not  be  important. 

It  is  the  purpose  now  to  develop  quantitative  criteria  for  neglecting 
either  one  or  the  other  effect. 

3,2  Deposition  from  polydispersed  cloud,  neglecting  diffusion 
The  presence  of  v  in  the  exponential  term  of  (15)  indi¬ 
cates  that  the  effect  of  varying  terminal  velocity  is  to  spread  the 
particles  much  like  a  diffusion  process  might  spread  particles  of 
uniform  terminal  velocities.  This  is  most  easily  seen  when  one 
considers  the  ground  deposition  distribution  pattern  resulting  from 
release  of  particles  whose  terminal  velocities  vary  like  (13),  in  a 
zero  diffusion  atmosphere.  For  this  case,  [neglecting  t]  in  (9)] 
a  particle  with  terminal  velocity  will  strike  the  ground  at 


(17) 


X  «  H/n 


1  • 

This  result  may  be  obtained  by  equating  the  right  hand  side  of  (4)  to 
zero,  and  letting  n  =  (2^  whence,  using  (6)  we  have 


hu(h)  H 

=  (i  +  a\n^  • 


This  expression,  in  which  the  distance  downwind  traveled  by  the 
particle  is  inversely  proportional  to  its  terminal  velocity  was  used  by 
Trouern-Trend  and  J.  Monaghan  (1956)  and  in  a  U.  S.  Discussion 
paper  presented  at  the  I2th  Tripartite  Conferenre  (1957). 

For  a  continuous  line  source,  the  deposit  per  per  sec  at 
2S  s  0  as  a  function  of  x,  resulting  from  the  lognormal  distribution  of 
mass  with  terminval  velocity  (13)  is 

(19*)  ■  exp 

Vzir  vx 

or 

(19b)  A^  =  ^^=^^exp^-ffn(L)-0]2/(2v2)].  . 

where  the  zero  subscript  indicates  deposition  for  zero  atmospheric 
diffusion  conditions.  If  the  terminal  velocity  of  the  particles  were 
constant,  all  the  particles  would  strike  the  ground  at  the  point  pre¬ 
dicted  by  (17).  The  effect  of  the  distribution  of  terminal  velocities 
is  to  spread  the  ground  deposition  over  a  wide  area.  The  effect  is  very 
much  like  a  diffusion  process,  even  though  in  this  case  atmospheric 
diffusion  has  been  assumed  to  be  zero.  The  distribution  predicted 
by  (19)  has  a  maximum  at 
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and  the  value  of  the  maximum  deposit  is  given  by 
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where  fl  is  defined  by  (14). 


or 


A  =-^=_L_e^4v' 


3.3  Analytic  approximation  for  polydispersed  clouds 

Since  the  varying  terminal  velocity  acts  as  a  diffusing  agent, 
it  seems  desirable  to  incorporate  this  effect  into  the  diffusion  equation. 
The  spreading  effect  of  the  varying  terminal  velocity  is  implicit  in 
the  term  w'c'  of  equation  (4).  If  the  diffusing  particles  have  widely 
varying  terminal  velocities,  u'  can  now  be  interpreted  as  the  devi¬ 
ation  of  the  velocity  of  the  particle  from  the  average  terminal  velo¬ 
city  of  all  the  particles  in  a  given  volume. 

Depending  on  the  relative  magnitude  of  v  and  p,  u>'  can  easily 
take  on  values  as  great  as  or  greater  than  w'.  For  this  case  we  con¬ 
sider  the  more  complete  diffusion  equation. 


We  now  identify  u'  as  a  characteristic  fluctuating  terminal  velocity 
which  is  proportional  to  the  arithmetic  standard  deviation  of  the  termi¬ 
nal  velocity  distribution.  Thus 

(23)  ^7?  =  kjcrz  _ 

where  k^  is  a  constant.  Thus  far  n*  in  equation  (22)  has  not  been 
defined  but  it  should  be  proportional  to  the  arithmetic  average  of  the 
terminal  velocities.  The  constant  of  proportionality  will  be  established 


below.  It  should  be  noted  at  this  po^t  that  since  the  heavier  particles 
fall  out  near  the  origin,  the  effective  H  and  o*  undoubtedly  decrease 
with  X}  it  is  not  certain  how  the  ratio  c/R  behaves  as  a  function  of 
distance  downwind.  In  any  event,  for  mathematical  feasibility  we  ne¬ 
glect  the  variation  of  R  and  a  with  x,  with  the  hope  that  adjustment 
of  the  two  constants  will  grossly  take  into  account  the  effect  of  this 
variation. 

By  using  (23)  in  (22)  it  can  be  seen  that  (22)  is  formally  the 
same  as  (5)  or 


= k  *  If + • 

The  only  difference  is  that  the  diffusion  coefficient  has  been  increased 
by  ah  amount  proportional  to  the  arithmetic  standard  deviation  of  the 
terminal  velocities.  The  solution  will  therefore  have  the  same  form 
as  (7)  and  may  be  written 


(24) 


with 


f*  _  H _  _  H 

(1  +  a)(ku»  +  kjcr)  "  T)  + 

_  n* 

^  “  (1  +  a)(ku*  +  k^tr)  “  r|  +,riQ  ’ 
where,  for  typographical  convenience,  we  have  let 
(25)  =  (1  +  a)kj(r 

Differentiating  (24)  and  changing  to  the  dimensionless  units  A’ 


and  x/f  of  (7),  we  have 
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It  remains  now  to  fix  the  constants  kj  and  0*.  This  is  done 
mo^t  simply  by  considering  the  solution  of  (22)  for  zero  diffusion 
conditions,  the  case  for  which  (19)  is  formally  valid.  For  u*=  0,  the 
solution  is  (24)  with  f)  =  0.  Thus 
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where  the  zero  subscript  indicates  a  zero  diffusion  solution.  In  particular 
the  point  of  maximum  concentration  and  the  value  of  the  maximum  concen¬ 
tration  for  r|  a  0  is 
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It  is  clear  that  since  the  variation  of  Rq  and  or  with  x  has  been  neglected, 
(28)  cannot  be  identical  with  (19),  but  since  we  have  two  undetermined 
constants  at  our  disposal  we  can  force  (28)  to  reproduce  the  location 


and  value  of  the  maximum  concentration  predicted  by  (19).  By  equating 


(29)  and  (20),  we  obtain 


”0  *  "o  "  =  S  (4^)’ 


and  by  equating  (30)  and  (21)  we  have 


ni  +  Po)L  «  J 


It  follows  from  (14),  (28).  and  (31)  that 


1  +  Pq  = 
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Using  (31)  in  (32)  ,  we  also  have 


Po/(i  +Po)  ri  +  Pof 

r(i  +  Po)  L  ®  J 


It  is  evident  that  pQ,  (1  +  a)kj  and  Oq/R  are  functions  of  v  only.  Equation 
(35)  is  relatively  easy  to  solve  for  1  +Pq  as  a  function  of  v.  The  solution 
is  graphed  in  Figure  1  as  a  function  of  v.  It  may  be  seen  that  p^  approaches 
zero  as  v  approaches  infinity  and  becomes  quite  large  as  v  becomes  small. 
Since  Pq  contains  the  diffusion  term  due  to  terminal  velocity  dispersion 
in  the  denominator,  it  is  clear  from  Figure  1  that  small  values  of  p^ 
imply  a  relatively  large  diffusion  constant  for  large  values  of  v  which  is 
a  quite  reasonable  result. 

Although  not  necessary  in  what  follows,  it  is  interesting  to  note 
the  behavior  of  (1  +  a)k^  and  kg  s  a*  functions  of  v.  This  is  shown 


Figure  2.  (1  +  a)k,  and  k_  as  functions  of 
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in  Figure  2  and  it  is  at  once  apparent  that  (1  +  a)k^i  which  may  be  re¬ 
garded  as  a  weighting  function  for  the  arithmetic  standard  deviation  of 
the  particle  terminal  velocity.becomes  large  as  v  becomes  large.  For 
small  V,  (1  +  a)kj  approaches  zero  and  approaches  1  indicating 

that  when  v  is  small  the  only  important  term  is  . 

The  constants  k^  and  k^  have  now  been  determined  from  the 
zero  diffusion  solution.  Thus  it  remains  to  incorporate  these  constants 
into  (24).  It  will  be  noted  that  is  still  largely  undefined.  If  we 
now  consider  the  case  where  v  is  very  small,  so  that  the  spread  of 
particle  terminal  velocities  may  be  neglected,  it  seems  reasonable  to 
assume  that 

4  —  — 

n  =  fi  and  p  =  . 

This  procedure  gives  us  the  correct  solution  when  there  is  zero  dis¬ 
persion  of  particle  terminal  velocity.  It  is  desirable  to  define  n* 
so  that  it  will  approach  SI  for  zero  terminal  velocity  dispersion  and 
will  approach  SIq  for  zero  atmospheric  conditions.  Such  an  expression 
is 


(36) 
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In. terms  of  the  statistics  of  the  terminal  velocity  distribution  and  arrang¬ 
ing  terms  so  that  the  variables,  ^  and  T|  of  (15)  and  (16)  are  introduced, 
we  have 
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and 
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3.4  Comparison  of  analytic  approximation  with  numerical 
integration  and  zero  diffusion  model 

The  solution  (24)  with  the  definitions  (36),  (37),  and  (38) 


constitute  an  analytical  solution  of  the  problem  we  posed  initially.  Several 
approximations  have  been  made  in  the  derivation,  and  it  is  desirable  to 
compare  the  predictions  of  (24)  with  the  results  of  the  numerical  inte¬ 
gration  of  (15)  which  is  certainly  more  exact.  In  fact  what  has  been 
done  thus  far  is  simply  to  invent  a  diffusion  process,  based  in  part  upon 
physical  notions,  which  it  is  hoped  will  reproduce  the  main  features  of 
the  numerical  integration.  Of  course  one  could  use,  for  n*  the  median, 
mode  or  some  other  measure  of  central  tendency.  This  would  give 
similar  but  not  identical  results.  ' 

Accordingly,  (15)  was  integrated  numerically  on  an  IBM  650 
computer.  The  limits,  zero  to  infinity,  were  replaced  by 
and  p  =  By  trial  and  error  it  was  found  that  dividing  p  into 

25  equally  spaced  intervals  was  siifficient  to  insure  adequate  accuracy 
in  the  numerical  results.  The  numerical  integrations  were  performed 
for  the  following  values  of  ^  and  v ; 
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The  numerical  integrations  are  shown  in  Figures  3  through  6  as 
indicated  above.  For  comparison,  the  analytic  approximation,  (34). 
with  the  coefficients  (36),  (37),  and  (38),  is  shown  along  with  the  zero 
diffusion  solution,  (19).  The  case  ^  s  •  1  (Fig.  3)  represents  the 
solution  when  the  atmospheric  diffusion  velocity  is  about  2.3  times 
greater  than  the  average  terminal  velocity  of  the  particles.  It  may  be 
seen  that  for  v  =  0.55,  the  solution  neglecting  atmospheric  diffusion 
overestimates  the  point  of  maximum  concentration  by  a  factor  of  about 
3,  and  the  value  of  the  maximum  concentration  by  a  factor  of  about  1/3. 

On  the  other  hand,  the  analytic  approximation  is  within  10%  of  the  lo¬ 
cation  and  value  of  the  maximum  concentration  and  fits  the  numerical  inte¬ 
gration  curve  reasonably  well  over  the  entire  range  of  integration.  It 
may  also  be  seen  that  for  the  case  v  a  1,  the  analytic  approximate 
underestimates  the  value  of  the  maximum  concentration  by  less  than 
10%  .  The  solution  neglecting  atmospheric  diffusion  yields  a  reasonable 
value  of  the  maximum  concentration  which  is,  however,  displaced  down¬ 
wind  by  more  than  a  factor  of  2.  Attention  is  also  called  to  the  fact  that 
the  analytic  approximation  is  very  close  to  the  numerical  integration 
past  the  point  of  maximum  concentration  where  the  solution  neglecting 
atmospheric  diffusion  is  particularly  bad. 

Figure  4  summarizes  the  computations  for  ^  =  0,  that  is 
the  case  where  the  atmospheric  diffusion  velocity  r)  is  equal  to  the 
median  terminal  velocity  of  the  diffusing  particles.  It  may  be  seen 
that  the  analytic  approximation  is  quite  satisfactory  for  all  three  values 
of  V.  The  solution  neglecting  atmospheric  diffusion  is  quite  bad  for 

V  =  0.30  and  0.55  and  yields  a  more  reasonable  approximation  when 

V  =  1. 


Figure  5  summarizes  the  results  for  (p  =  0.81  and  1.50.  Again 
the  analytic  approximation  is  quite  close  to  the  numerical  integration 
while  the  solution,  neglecting  atmospheric  diffusion,  although  closer 
to  the  numerical  integration  than  was  the  case  for  small  values  of  0, 
is  still  not  as  good  as  the  analytic  approximation. 

Finally,  Figure  6  summarizes  the  results  for  p  =  1.6i  and  2.30. 

It  will  be  noted  that  for  these  large  values  of  p  (i.  e.  ,  equals  about  5 
or  10  times  r|),  the  solution  neglecting  atmospheric  diffusion  is  toler¬ 
able  for  v  =  0.53.  The  analytic  solution,  however,  is  still  a  better  ap- 
proximation  even  for  these  large  relative  values  of  H. 

In  summary,  the  analytic  approximation  fits  the  more  exact 
numerical  integration  reasonably  well  for  all  values  of  p  and  v  inves¬ 
tigated.  It  may  be  seen,  however,  that  the  solution  neglecting  atmo¬ 
spheric  diffusion  is  quite  bad  for  small  values  of  (p  and  v  but  becomes 
tolerable  as  p  and  v  become  large.  This  is  not  a  surprising  result 
since  a  large  value  of  <p  and  v  imply  that  atmospheric  diffusion  process¬ 
es  are  unimportant  compared  to  the  median  terminal  velocities  and  the 
spread  of  terminal  velocities  around  the  mean. 

4.  Relative  importance  of  diffusion  velocity  and  the  spread  of  terminal 
velocity 

4. 1  When  can  atmospheric  diffusion  be  neglected? 

Now  that  the  analytic  approximation  to  the  integral  has  been 
established  to  be  reasonably  accurate,  it  is  possible  to  develop  simple 
expressions  for  estimating  under  which  circumstances  atmospheric  dif¬ 
fusion  and/or  variations  in  particle  terminal  velocity  may  be  neglected. 
For  example,  equation  (21)  represents  a  solution  in  which  atmospheric 


diffusion  is  completely  neglected.  Under  what  circumstances  is  this 
neglect  justified?  The  question  is  readily  answered  if  we  define  a 
quantity  R*  which  is  the  ratio  of  the  location  of  the  maximum  concen¬ 
tration  predicted  by  equations  (20)  or  (29)  which  neglects  atmospheric 
diffusion  to  that  predicted  by  (26)  which  includes  the  effect  of  both 
atmospheric  diffusion  and  dispersion  of  particle  terminal  velocities. 
Similarly  a  quantity  A*  relating  the  ratio  of  the  predicted  maximum 
deposit  for  the  zero  diffusion  assumptions  (21)  or  (30)  to  that  pre¬ 
dicted  by  (27).  Thus 
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Equation  (39)  is  graphed  in  Figure  7.  If  R*  is  close  to 
unity  there  is  little  difference  between  the  location  of  the  maximum 
concentration  predicted  by  either  equation.  It  may  be  seen  that  neglect 
of  the  atmospheric  diffusion  term  generally  leads  to  an  overestimate 
of  the  location  of  the  maximum  concentration.  It  may  also  be  seen 
that  for  any  value  of  p,  there  always  exists  a  v  for  which  atmo¬ 
spheric  diffusion  may  be  neglected  with  tolerable  results.r,  For  ex¬ 
ample,  the  error  in  the  location  of  the  maximum  concentration  caused 
by  neglect  of  atmospheric  diffusion  is  10%  or  less  for  ^  >  2  and 
V  >  0.9 ;  for  ^  =  0,  i.  e.  ,  t)  =  W  ,  the  error  is  about  80%  for  v  =  0.55, 
about  45%  for  v  s  1,0  and  about  127o  for  v  =  1.5  .  Since  the  coefficient 
of  variation  of  the  distribution  of  mass  as  a  function  of  terminal  velocity 


is  (e'^  -  1)  ,  a  value  of  1.5  for  v  implies  that  (r/il  *  2.9  which 

represents  a  very  wide  distribution  of  terminal  velocity. 

The  curve  of  A*  as  a  function  of  0  and  v  is  shown  in  Figure  8. 
Again  the  zero  diffusion  solution  badly  overestimates  the  value  of  the 
maximum  concentration  for  small  values  of  v.  The  percent 
error,  however,  becomes  quite  small  (less  than  10%)  for  v  =  0.65, 

0  =  -1.0;  V  a  0.77  for  0  =  0,  etc.  In  evaluating  these  numbers  it 
should  be  borne  in  mind  that  we  are  dealing  here  with  an  infinite  line 
source  solution,  and  in  general  all  errors  in  concentration  would  be 
roughly  squared  when  applied  to  a  point  source. 

From  inspection  of  (39)  and  (40)  it  is  apparent  that  since 
Pq-*P*-*’  0  as  V -•  00  ,  both  R*  and  A**^l  for  large  values  of  v. 

The  rate  at  which  both  ratios  approach  one  differs.  For  example, 
at  V  =  1,  0  =  -1,  the  estimate  of  the  point  of  maximum  concentration 
is  off  by  a  factor  of  2,  but  the  value  of  the  maximum  concentration  is 
approximately  correct.  This  case  was  exhibited  in  Figure  3.  It  seems 
clear  that  the  solution  neglecting  atmospheric  diffusion  predicts  the 
value  of  the  maximum  concentration  for  relatively  large  v  much  better 
than  it  predicts  the  location  of  the  maximum.  Typically,  this  yields 
curves  of  deposition  which  peak  further  downwind  than  the  correct 
solution,  and  rather  substantial  differences  in  predicted  deposition 
result- -sometimes  for  considerable  distances  past  the  point  of  maxi¬ 
mum  concentration.  This  occurs  even  when  (30)  predicts  the  value  of 
the  maximum  concentration  correctly.  For  this  reason  is  is  best  to 
use  R*  or  Figure  7  as  the  criteria  rather  than  A* 


(46) 


Inspection  of  (43)  and  (44)  indicates  that  everything  else  being 
held  constant,  the  percent  error  in  the  location  of  the  maximum  concen¬ 
tration  is  less  when  the  mean  rather  than  median  terminal  velocity  is 
employed  in  the  definition  of  p.  Equation  (43)  is  graphed  in  Figure  9  as 
a  function  of  v  and  It  maybe  seen  that  the  error  in  the  location  of 
the  maximum  ground  deposit  is  negligible  when  v  is  Small '!• -as  indeed  it 
should  be.  The  percent  error  becomes  large  as  v  increases,  but  the 
results  are  still  tolerable  (less  than  15%  for  v  <  0.55  and  p  <  0.5)  for 
surprisingly  large  values  of  v  . 

Equation  (45)  is  sketched  in  Figure  10.  The  value  of  the 
maximum  concentration  is  apparently  more  sensitive  to  the  spread  of 
terminal  velocity  than  is  the  point  of  maximum  concentration.  For 
fixed  values  of  terminal  velocity,  and  for  small  values  of  diffusion  velo¬ 
city  (0  positive  and  large)  it  is  apparent  that  the  spread  in  terminal 
velocity  cannot  be  neglected  even  when  v  is  as  small  as  0.3.  However, 
when  the  diffusion  velocity  is  large  compared  to  the  mean  terminal  velo¬ 
city  (p  negative),  the  dispersion  of  terminal  velocity  may  be  safely 
neglected.  Thus  when  b  =  -1,  the  error  in  and  x  is  less  than 
15%  for  V  as  large  as  0.8. 


5.  Atmospheric  Applications 
5. 1  Neutral  conditions 

It  has  been  shown  that  the  behavior  of  the  ground  deposition 
resulting  from  the  release  of  particles  whose  mass  distribution  is  log> 


Figure  9-  as  a  function  of  v  for  indicated  values  of  ^  . 
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normal  with  respect  to  terminal  velocity  depends  upon  the  values  of  two 
parameters.  The  first  parameter  p  represents  the  logarithm  of  the 
ratio  of  the  mass  median  terminal  velocity  to  the  virtual  diffusion  velo¬ 
city,  T).  The  second  parameter,  v  .  depends  on  the  dispersion  of  the 
terminal  velocities  around  the  logarithmic  mean.  The  only  meteorolo¬ 
gical  pars^meter  is  T)  =  (1  +  a)ku*.  Here  (1  +'a)  is  a  measure  of  the 
wind  shear  which  varies  with  atmospheric  stability  and  ground  rough¬ 
ness.  The  quantity  u*  is,  of  course,  proportional  to  the  square  root 
of  the  surface  stress  which,  in  turn,  is  a  function  of  wind  speed,  ground 
roughness  and  stability. 

The  solutions  presented  thus  far  are  valid  for  an  adiabatic 
atmosphere.  If  (1  +  a)  is  defined  to  yield  the  same  average  wind  speed 
between  ground  level  and  height  h  and  the  same  speed  at  release  height 
h  as  that  predicted  by  the  more  accurate  logarithmic  law,  it  follows 
that 


(47) 


(1  +  a)  = 


fn(h/z^) 
<n(h/ z^)  - 1 


(48) 


and 

(■49) 


n[fn(h/zj-lj 

P  =  fn  - % - 

u(h)k 


The  parameter  p  is  therefore  a  function  of  (h/z^),  P  and  u{h).  The 
solid  line  plots  in  Figures  11  and  12  are  curves  of  constant  p  as  functions  of  u(h) 
and  £2  for  values  of  h/z^  =  1.5  x  10  and  1.5  x  10  respectively.  In  terms 
of  our  previous  work  these  values  of  h/z^  correspond  to  a  release  height 
of  15  m  and  a  z^  of  1  cm  and  10  cm  respectively. 


It  may  be  eeen  that  for  fixed  wind  speed  and  mass  median 

terminal  velocityi  the  value  of  p  decreases  as  increases.  If  w 

—  2 
is  60 cm/sec  and  u(h)  is  5  m/sec,  ^  1.2  and  1.7  for  h/ z^  a  1.5x10 

3 

and  1.5  x  10  respectively.  Referring  to  Figures  7  and  8,  it  may  be 
seen  that  for  these  values  of  p  ,  atmospheric  diffusion  may  be  ne¬ 
glected  with  fair  approximation  for  values  of  v  >  0.6  for  the  smooth 
terrain  and  0.7  for  the  rougher  terrain.  Similarly,  it  may  be  seen 
from  Figures  9  and  lO,  that  for  these  values  of  u  and  11  ,  the  effect 
of  dispersion  of  terminal  velocity  on  the  diffusion  pattern  cannot  be 
neglected  except  for  values  of  v  less  than  03.  If  ^  remains  constant 
and  the  wind  speed  at  release  height  increases  to  10  m/sec,  ^  ranges 
between  0.4  and  0.9  and  the  range  of  validity  of  the  two  approximations, 
one  neglecting  atmospheric  diffusion  and  the  Other  neglecting  terminal 
velocity  dispersion,  will  change  accordingly.  It  is  interesting  to  note 
that  the  case  12  =  30  cm/sec,  u  =  5  ih/sec  corresponds  to  the  same 
as  f2  a  60  cm/sec,  u  =  10  m/sec. 

This  latter  point  illustrates  the  efficiency  implicit  in  the 
non-dimensional  plots  of  Figures  3-6.  From  the  definition  of  (48),  it 
is  evident  that  a  large  number  of  possible  12,  u(h)  and  h/z^  combi¬ 
nations  can  combine  to  yield  the  same  value  of  p.  For  any  such  com¬ 
bination  yielding  a  fixed  value  of  say  ^  =  1.61  and  a  v  value  of 
0.53,  the  appropriate  curve  of  Figure  6  would  be  valid.  The  dimensional 
deposition  may  be  recovered  by  multiplying  the  ordinate  by  Q/f.  For 
an  adiabatic  atmosphere,  using  the  definition  of  (1  +  a)  introduced  at 
the  beginning  of  this  section,  it  follows  that 


'  h.[in(h/zj- 1]^ 

(50)  ^  =  Til  '  Z,  ,  ;  .  *  gm/(m seq) 

If,  for  example,  the  height  of  release  were  15  m,  and  z^sO.O'lm, 
then  f  would  have  the  value  of  511  m,  and  the  dimensional  deposition  in 
gni/(m^sec)  for,  say,  Q=1000 gm/(m sec)  would  be  recovered  by  multiplying 
the  ordinate  of  the  appropriate  figure  by  Q/f  =  1.96.  The  dimensional 
distance  may  be  recovered  by  multiplying  the  x/f  scale  by  511  m.  The 
subsidiary  scales  on  Figure  6  indicate  the  values  thus  obtained.  This 
case  corresponds  closely  to  the  example  presented  in  Progress  Repor.t 
523.19  for  £1  =  e^'®^  =  0.58  m/aec  and  u(h)  =  4.6  m/sec. 

5.2  Diabatic  conditions 

An  approximate  method  of  extending  the  monodispersed 
particle  solution  (7)  to  the  diabatic  atmosphere  has  been  proposed  by 
Godson  (1958).  For  the  diabatic  atmosphere  it  is  assumed,  following 
Deacon  (1949).  that 


where  8  varies  with  stability  in  such  a  fashion  that 
8  <  1,  for  stable  conditions 
8  =  1,  for  adiabatic  conditions 
8  >  unstable  conditions  . 

It  has  long  ago  been  observed  (Davidson  and  Barad,  1956;  Panofsky, 
Blackadar  and  McVehil,  I960;  and  others)  that  the  8  formulation  is 
qualitatively  valid,  but  quantitatively  incorrect.  In  view  of  the  approxi- 


mations  that  follow,  the  whole  subsequent  argument  is  seml-qu^ntitdtive 
in  nature,  and  probably  does  not  warrant  dealing  with  the  somewhat 
more  refined  but  nevertheless  incorrect  (for  stable  conditions)  versions 
of  the  log -linear  law.  or  the  Ellison  interpolation  formula  contained  in 
the  ^ast  referenced  paper. 

At  any  rate,  (5)  and  (22)  cannot  be  solved  in  terms  of  tabulated 

a 

functions  unless  K  is  a  linear  function  of  height.  Since  is  linear 
only  for  p  =  1  (adiabatic  conditions).  Godson  suggested  that  an  ap¬ 
proximate  way  to  deal  with  the  problem  is  to  find  an  t  such  that 


It  follows  that  when 


«  -  2  /Vj-P 


the  mean  value  of  the  diffusion  coefficient  from  zero  to  h  under  any 
stability  conditions  is  reproduced  by  the  linear  diffusion  coefficient 

(53)  K  =  ku*<z  . 

Accordingly  an  approximate  solution  of  (5)  or  (22)  for  diabatic  con¬ 
ditions  may  be  written  exactly  in  the  form  (7)  or  (24)  provided 

(54)  tip  =  (1  +  o)kfu’'‘  , 


From  (52),  it  is  apparent  that 
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and 

where 

-  1 

(55)  U=^/u(h)dh 

The  parameter  (p  is  now  determined  from  the  above  as  a  function  of 
h/z^,  p  and  u(h).  The  dashed  lines  on  Figures  11  and  12  are  lines  of 
constant  <p  for  p  =  1.10.  Experience  indicates  that  this  value  of  p  cor¬ 
responds  to  a  case  of  marked  instability  and  it  is  unlikely  that  the  wind 
speed  at  release  height  will  be  greater  than  5  m/sec  for  this  value  of  p. 

It  is  apparent  that  the  unstable  atmosphere  decreases  the  value  of  p 
from  its  adiabatic  value.  For  example,  for  h/z^  =  1.5  x  10  , 

=  60  cm/sec,  u(h)  =  5  m/sec,  <p  decreases  from  about  1.2  to  about 
0.5.  Consequently,  neglect  of  atmospheric  diffusion  is  justified  only  for 
larger  values  of  v  than  was  the  case  for  the  adiabatic  atmosphere. 
Similarly,  neglect  of  the  effect  of  the  variation  of  terminal  velocity  on 
the  ground  deposition  pattern  would  be  justified  for  slightly  larger  values 
of  V  than  was  the  case  for  the  adiabatic  atmosphere. 

The  dotted  lines  on  Figures  11  and  12  represent  constant  values 
of  0  for  p  a  0.8.  This  value  of  p  corresponds  to  a  moderately  strong 
inversion.  It  may  be  seen  that  the  values  of  0  are  increased  considerably 
over  the  adiabatic  values.  Consequently,  during  inversion  conditipns, 


(1-  P)ku(h) 

-  i 


the  primary  diffusing  agent  is  the  variation  in  terminal  velocity,  and 
neglect  of  atmospheric  diffusion  is  valid  foi;  relatively  small  values  of 
V  .  On  the  other  hand,  a  model  which  employs  atmospheric  diffusion 
and  the  average  terminal  velocity  only,  would  be  quite  bad  for  even 
small  values  of  v  . 
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